1/ 


AD- A 1 46  265  THE  LINDSTROM  MADDEN  METHOD  FOR  SERIES  SYSTEMS  WITH 
REPEATED  COMPONENTS(U)  WISCONSIN  UNIVMAD1SON 
MATHEMATICS  RESEARCH  CENTER  A  P  SOMS  AUG  84 
UNCLASSIFIED  MRC-TSR-2732  DAAG29-80C-0041  F/G  12/1  NL 


MICROCOPY  RESOLUTION  TEST  CHART 
naiional  eunuu  01  atandarcs  rjt j  a 


! 

i 


Mathematics  Research  Center 
University  of  Wisconsin— Madison 
610  Walnut  Street 
Madison,  Wisconsin  53705 


August  1984 

(Revived  July  25,  1984) 


Sponsored  by 

U.  S.  Army  Research  Office 
P.  O.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 


Approved  (or  public  release 
Distribution  unlimited 


Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  VA  22217 


UNIVERSITY  OF  WI9C0NSIN-MADIS0N 
MATHEMATICS  RESEARCH  CENTER 

THE  LINDSTROM-MADDEN  METHOD  FOR  SERIES  SYSTEMS  WITH  REPEATED  COMPONENTS 

Andrew  P.  Some* 

Technical  Summary  Report  #2732 
August  1984 

ABSTRACT 

The  Lindstrom-Madden  method  of  computing  lower  confidence  limits  for 
series  systems  with  unlike  components  is  extended  to  series  systems  with 
repeated  components  utilizing  the  results  of  Harris  and  Soms  (1983).  An  exact 
solution  is  given  for  no  failures  and  key  test  results,  together  with  an 
approximation  for  the  general  case.  Numerical  examples  are  also  provided. 
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SIGNIFICANCE  AND  EXPLANATION 


Series  systems  with  repeated  components  arise  often  in  engineering  and 
physics.  It  is  therefore  important  to  utilize  data  obtained  on  individual 
components  in  an  efficient  manner  when  assessing  the  reliability  of  the 
combined  system.  This  paper  gives  o  .u  method  for  doing  so. 


The  responsibi lity  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MR C,  and  not  with  the  author  of  this  report. 


THE  LIND3TR0M-MADDEN  METHOD  FOR  SERIES  SYSTEMS  WITH  REPEATED  COMPONENTS 


Andrew  p.  Sons* 

1.  INTRODUCTION  AND  SUMMARY 

A  problem  of  substantial  Importance  to  practitioners  in  reliability  is  the 
statistical  estimation  of  the  reliability  of  a  series  system  of  stochastically  independent 
components  when  some  components  are  repeated,  using  experimental  data  collected  on  the 
individual  components.  In  the  situations  discussed  in  this  paper,  the  component  data 
consist  of  a  sequence  of  Bernoulli  trials.  Thus,  for  component  i,  i  ■  1,2, ...,k,  the 
data  is  the  pair  (n^Y^,  where  n^  is  the  number  of  trials  and  Ya  is  the  number  of 
observations  for  which  the  component  functions.  Yj  ,Y,,  •  •  •  ,Yj.  are  assumed  to  be  mutually 
independent  random  variables.  We  assume  that  there  are  components  of  type  i, 

1  <  i  <  k.  Then  the  parameter  of  interest  is  h (p^ ,pj, . . . ,p^l  ■  h(p),  the  reliability  of 
the  system,  whei  e 

k  f 

h(p)  -  I  I  p. 1  . 
i-l  1 

More  specifically,  it  is  desired  to  obtain  a  Buehler  (1957)  optimal  lower  1  -  a 
confidence  limit  on  h(p). 

The  case  of  Y1  -  y2  «  . . .  «  -  1  has  been  treated  in  Sudakov  (1974),  Winterbottom 

(1974),  and  Harris  and  Soms  (1983). 

In  Section  2  we  summarize  the  general  theory  of  Harris  and  Soms  (1983)  applicable 
here.  In  Section  3  the  exact  solutions  to  no  failures  and  key  test  results  are  given. 
Lindatrom-Madden  type  approximations  are  given  in  Section  4.  Section  5  contains  numerical 
examples. 
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2.  BUEHLER'S  METHOD  FOR  OPTIMAL  CONFIDENCE  LIMITS 
We  now  specialize  the  general  results  of  Harris  and  Soms  (1983)  on  optimal  confidence 
limits  for  system  reliability  to  a  series  system  with  independent  and  repeated 
conqponents .  As  in  Section  1,  let 

„  k  y. 

h(p)  =■  TT  p.  , 

i=1 

o  <  <  1,  -  n±  -  Yj,  «  nt  -  yif  1  <  i  «  k,  S  -  (x^  *  0,1,...,^,  1  <  i  <  k}  and 

let  g(x)  =  (x1,x2»...,x){)  be  an  ordering  function,  i.e.,  for  real  xit  0  <  xi  <  ni( 
g(x)  is  non-decreasing  in  each  component.  It  is  often  convenient  to  normalize  g(x)  by 
letting  g(0)  =  1  and  g(n)  =  0.  With  such  a  normalization,  g(x)  is  often  selected  to 
be  a  point  estimator  of  h(p).  Also  let  R  =  ^rl'r2'”*'r  '  s  >  2}  be  the  range  set  of 
g(x).  with  no  loss  of  generality  we  order  R  so  that  r-]  >  r2  >  ...  >  rs  and  let 
“  {x|g(x)  =  ri(  X  e  S,  i  =  1,2,.. .,s}.  The  sets  A^  constitute  a  partition  of  S 
induced  by  g(x).  We  assume  throughout  that  the  data  is  distributed  by 


f (x;p) 


PMX 

P 


x) 


k  n.  y.  n.-y. 

FT  t/KV  1  . 

i«1  yi 


(2.1) 


where  g^  =  1  -  p^,  i  =  1,2,...,k.  with  no  loss  of  generality,  we  assume 
<  n2  <  ...  <  n^. 

From  these  definitions,  it  follows  that 

j 

pplx  '  U  =  P^{g(X)  >  r..}  .  (2.2) 


From  (2.1)  and  (2.2),  we  have 
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(2.3) 


P'{g(X)  >  r  }  =  i  l  •••  l  f(i;p) 

4r°  i2=°  v° 


where  i  »  (i  1 ,  i2 ,  • . . ,  i^ )  and  U2  =  u2(  if ) , . . .  ,uk  =  u^li  1  ,  i2 ,  . .  • ,  ijt_  ^ )  are  integers 
determined  by  ry  Equivalently, 


(t1]  tt2]  [tk] 


P' 

P 


Jg(X)  >  r.)  -  l  l  “•  l  f ( i ; p ) 
J  i  •;  =n  4  =r\ 


(2.4) 


V°  i2=o  ik=o 


where  t2  =  t2(i^  ), . . .  .tjj  =  tk ( i j ,i2, . . . , ik_ 2 ) ,  with  t^  =  sup{t!0  <  t  <  n-j  and 
g(t ,0 ,0 , . , . ,0 )  >  r^}  and  t  j(i.,  ,i2 , . . .  ,ij_  1 )  =  sup{t|o  <  t  <  n^  and 
g(i1,i2,...,i4_1,t,0,...,0)  >  r j ) >  t  =  2,3,. ...k. 

We  now  introduce  the  notion  of  Buehler  optimal  confidence  limits.  Let  g(x)  “  r^, 
Then  define 

ag(Jj  *  inf {hCp > |P^{i|g(i)  »  g(x)}  >  a}  .  (2.5) 

Equivalently,  by  (2.2),  we  can  also  write 


g(x) 


inf |h(p) |P'|x  e  U  A . }  >  o[  . 
P  1-4  * 


i=  1 


1  1 


Then  we  have,  from  Harris  and  Soms  (1983), 

Theorem  2.1.  a  is  a  1  -  a  lower  confidence  limit  for  h(p).  if  b  is  anv 

-  g( x )  c  g(X)  1 

other  1  -  a  lower  conf 


idence  limit  for  h(p)  with  b  >  b  >  ...  >  b  ,  then 

1  r2  3 


b  <  a  for  all  x  e  s. 
g(x)  g(x) 

Two  possible  choices  of  g(x)  are 


'  *  T  ■ 

g(x)  »  1  r  ((n4  -  x^/nj  1  , 

i“1 


(2.7) 


or 
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g<5) 


(2.8) 


Both  reduce  to  the  generally  used  g(x)  for  series  systems  with  independent  components 
when  >i  ■  *2  *  •••  "  1|5  “  1*  i.e., 

k 

g(x!  =  "1  f  (n^  -  x^)/"*  . 
i=  1 


Since  (2.7)  is  the  maximum  likelihood  estimator  of  h(p)  we  will  use  it  here  and  from  now 
on  it  will  be  understood  that  g(x)  is  given  by  (2.7).  With  this  choice  of  g(x),  we 
assume  from  now  on  that  0  <  x^  <  n^,  i  =*  1,2,...,k,  since  ag(J)  “  0  if  some 
xi  “  ni*  With  this  assumption,  the  t^  in  (2.4)  are  given  by 


and 


t1  “  n1 


V7, 


(2.9) 


n,  -  ITT  «n 

'  i-1 


t-1 

IT  (ns 

S*1 


(2.10) 


l  “  2,...,k,  with  1  T  n.  »  1. 

i=k+1  1 

For  the  purpose  of  simplifying  the  calculation  of  ag($)  in  special  cases  it  is 

necessary  to  state  additional  results  from  Harris  and  Soms  (1983). 

Theorem  2,2.  Let  g(x)  =  r^  and  let 

f*(xia)  =  sup  P'{g(X)  >  r  },  0  <  a  <  1  .  (2.11) 

h(p)=a  P  3 


Then 
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inf  f  (xia)  -  0,  sup  f  (x»a)  »  1 
0<a<  1  0<a<  1 


and  f  (xia)  is  strictly  increasing  in  a. 

*  „ 

Theorem  2.3.  f  (xja)  «  a  has  exactly  one  solution  a_  in  a  and  a  «  a_,^,. 
-  ■*  a  a  gix) 


3.  EXACT  SOLUTIONS  FOR  ZERO  FAILURES  AND  KEY  TEST  RESULTS 
We  first  assume  that  x  »  (0,0,.  ..,0)  •  0  and  use  Theorem  2.3  to  obtain  a^~j. 
Theorem  3.1.  If  x  <»  0,  then 


•  '  -rr  ni 

f  (Oja)  *  sup  ]  I  Pi 


K  yi  im1 


where  n./y .  «  min  n, /y .  and 
3  3  1<i<k  1  1 


a  ^  ■  a 

g(0) 


FTP?  -  ITT  Pj)  j  j  TT 


since  n^Y^  ~  n^Y^  *  0  is  equivalent  to  n^/y^  *  n^/Yj,  which  is  true,  and  therefore 

*  (niWi)A1  1Ai 

|lp,  J  J  J  <  1.  (3.1)  follows  by  noting  that  the  choice  p.=  »  a  ,  p,  »  1, 

i-1 


TT-  "i 

i  *  j,  gives  |  )  p,  “  a  .  Then,  using  Theorem  2.3,  we  obtain  (3.2),  which  reduces 


to  the  known  series  result  if  Y1*Y2“»**’*Yk“1. 


We  now  turn  to  analogues  of  key  test  results  (see,  e.g.,  Winterbottom  (1974)  and 

Harris  and  Soms  (1983)).  We  define  a  key  test  result  if  y.  =  max  y.  (recall  that 

14i4k  1 

n1  =  min  n^)  and  x  =  ( x.j  ,0 , . , .  ,0  ) . 

1<i<k 

Theorem  3.2.  If  x  is  a  key  test  result  and 


k  Y  k  Y  k 

1*|TT  <ni  -  1  >  TT  ("i  '  *L)  M  =  {z|  l  In^z^ 


i=1 


i=1 


i=*1 


>  y  (n. -x.  )} 

.  .  i  i  ' 

i=1 


(3.3) 


then 


f  ( x ;  a ) 


I  1/Y  (n  -  x,,*,  +  1) 
a 


(3.4) 


where  I„(a,b)  is  the  incomplete  beta  function.  Let  b  denote  the  solution  in  b  of 
A  a 

a  =  Ib(ni  -  x1(x,  +  D  . 

Then  ag(J)  =  ba  •  Note  that  bQ  is  the  usual  1  -  a  lower  confidence  limit  on  p, 
given  xi  failures  in  n-j  trials. 

Proof.  Without  loss  of  generality  we  can  assume  that  n1  =  n2  =  ...  =  nk,  for  otherwise 
we  can  write  (2.4)  as 


P-{g(X) 

P 


a  1  ?  ,n1,  "rS  S  Xl"  1  ,n2,  X2 

*  rj}  =  *  „  li JP1  q1  .  \  Jp2 

i  =0  1  i  =0  2 


Xr1ri2'-”_ik-2  n 


Vi-0 


/  k-1,  nk-l"Vl  k-1 

^k-l  k_1  qk_1  Ipk  nk" 


*  *  •  ” ^ k“  1  ^  > *  1  ^ t ) 


x.  x-i ,-i. . .  .-i,  _  .  . 

V1  ,n1  )  nrl1  X1  1  1  y  k"2fnk-1  ,  nk-1-ik-1  ik-1 

<  l  (4  )P,  9,  l  n  (,  JPfc-1  qk-1  V  ("l  ■ 


ir°  1 


ik_i=°  *k-i 


(x^*i^“i2_. • .-i^.i ) ,  x  i ^*i2> . )  , 


(3.5) 
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where  g(x)  -  rj,  by  the  monotone  likelihood  ratio  property  of  the  beta  distribution 
(Ix(a,b)  has  a  monotone  likelihood  ratio  in  -a  for  fixed  b,  which  implies  that 
Ix(a,b)  is  a  decreasing  function  of  a).  A  similar  argument  applies  to  the  other 
indexes.  Thus,  if  (3.4)  is  true  for  nj  »  n2  =  ...  -  n^,  by  (3.5)  it  follows  for 
n^  <  n,  <  ...  <  ny. 

So,  asstiming  n  «  (n1  ,n.  , . . .  ,n1 ) ,  we  seek  to  maximize 
k  n1  k  k 

l  1  >  1  (ni  -  xi>  ■  i  -  (3-6 

p  i-1  j-1  3  i-1  i-1 


where  .  are  independent  Bernoulli  random  variables  with  parameter  and 

k  y  k  Yi  k  1/Yj 

I  |  p,  «  a.  If  |  |  p.  »  a,  then  "|  T  Pi  ranges  from  a  to 

i-t  1  i-1  i-1 


,/Y1 

a  ,  r,  -  min  Y, •  This  is  seen  as  followst 
3  1<i<k 


k  k  Y,  1/Y .  k  1-Y,/Y.. 

TT  Pi -ITT  TT  Pi 

x-1  i-1  i-2 


1/Y,  k  (Y^Y.I/Y,  1/Y, 

a  TT  P,  <  « 

i-2 


and 


k  k  Y,  1/Yj  k  1-Y,/Y. 

TT  Pi  -  (TT  pt  )  j  TT  Pi 

i-1  i-1  i-1 


1/Y.  k  ( Y . _Y. )/Yj  1/Y4 

.  j  TT  Pi  3  j  »  «  3 

i-1 

i*1 


1/T, 


and  the  choices  p^  “  a 


P2 


1/Y. 


-  1,  and  Pj  »  a  ,  p,  ■  1,  i  1  j,  attain 


these  values.  From  the  results  of  pledger  and  Proschan  (1971),  for  each  b  -  1  f  pi( 

i-1 
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1/Y.  1/Y1 

a  2  <  b  <  a  ,  (3. 6)  is  maximized  by  p^  =  b,  p^  =  1 ,  2  <  i  4  k.  Further,  the 

Vy.  i/y,  Vy  , 

maximum  over  b,  a  4  b  <  a  ,  of  the  maxima  for  each  b  is  given  by  p,  =  a  , 
Pj  “  1,  2  <  i  <  k,  by  the  monotone  likelihood  ratio  property  of  the  binomial 


VY, 


distribution,  and  p1  =  a 
the  proof. 


'  Pi 


1,  2  <  i  <  k,  satisfies 


k  Yt 

p^  =  a.  This  completes 


i=1 


If  i  =  y2  =  • • •  =  ^k  =  1 '  some  guidelines  for  the  verification  of  (3.3)  are  given 

in  Harris  and  Soras  (1983).  In  the  present  case  (3.3)  must  be  verified  by  trial  and  error 

k  Y .  Y .  k 


by  showing  that  min 

k  i=1 

I  x,»x 


k  Y i  Y,  k  Y 

TT  (ni  -  x  )  =  (n,  -  x1  )  TT  nA 

i=1  i=2 


and  that 


i=1 


1 


k  Y .  Y1  k  Y  ■ 

max  "j  T  ("i  -  x  )  <  ("i  *  X  )  1  T  n  • 

k  i=1  i=2 

)  x.»x,+1 

i-1  1  1 

Example  3.1.  Let  k  =  3,  n  =  (5,5,5),  Y  =  (3,3,2),  a=  .10  and  x  =  (1 
3  Yi  3  Y± 

min  "|  J"  (nA  -  x^)  1  =  200000  and  max  TT  (nA  -  x  >  =  1406. 

3  i=1  1  3  i=1 

l  *i=1  1  *i=2 

i=*1  i=1 

x  is  a  key  test  result  and  (3.3)  is  satisfied  and  hence 


) .  Then 


where 


.10 


a  »  .4161-5  »  .0720  , 
g(x ) 

I.4161(4'2'’  Further,  it  can  also  be  verified  that  x  =  (2,0,0)  is  a  key 


test  result  for  which  (3.3)  is  satisfied,  but  that  for  x  =  (3,0,0),  (3.3)  is  violated. 

7, 


Note  that  Theorem  3.2  asserts  that  a  =  b 

g(x)  a 


for  0  <  a  <  1.  It  is  thus  possible 
that  (3.3)  is  not  true  but  the  conclusion  still  holds  for  a  of  practical  importance. 

This  is  taken  up  in  Section  4. 


4.  THE  LINDSTROM-MADDEN  METHOD  FOR  SERIES  SYSTEMS  WITH 
REPEATED  COMPONENTS 

When  Y1=Y2='-*=Yr=V  the  Lindstrom-Madden  method  (henceforth  abbreviated 


L-M)  is  an  approximation  b^-^  to  ag(x)  t*le  ^orm 


bg(x) 


min 

14i<k 


b  (n  . )  , 
a  l 


(4.1) 
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where 


a 


Zb  <n  )<ni 
a  i 


t0i't0i  +  **  ' 


(4.2) 


with 


t0i  «  nt(1  -  T~T  ("i  -  xi)/ni)  , 
i«1 


(4.3) 


i.e.,  tQ1  is  the  maximum  of  the  recursive  indexes  t^  defined  by  (2.4).  For  the  usual 
levels  of  a,  ( x )  "  ba*nl'"  Purther,  numerical  evidence  indicates  (Harris  and  Soms 
(1983))  that  for  a  levels  of  practical  significance 

(4.4) 

(4.4)  was  incorrectly  claimed  to  be  true  for  0  <  a  <  1  in  Sudakov  (1974)  and  this  is 
discussed  at  length  in  Harris  and  Soms  (1983).  However,  (4.4)  is  known  to  hold  for 
special  cases  (Winterbottom  (1974)  and  Harris  and  Soms  (1983)). 


bg(x)  <  ag(x) 


Motivated  by  the  above,  we  now  give  an  L-M  approximation  b  .*•,  to  a  ~  for 

9'X J  9' ^  ' 

arbitrary  by 


bg(x) 


min  b  (n  .  ) 
1<i<k  °  1 


(4.5) 


where 


a  =  xb  (n,  )(ni  “  t0i't0i  +  1 '  ’ 
a  i 


(4.6) 


with 


_!<_  Y,  /  k  Y,  VY. 

fc0i  *  ni  ”  U  I  ~  xi>  /  TT  n  J  ,  (4.7) 

j“1  3  /  j“1  3 

j*i 

i.e.,  tg^  is  the  maximum  of  the  recursive  indexes  tj  defined  by  (2.4).  However,  in 
this  case  it  is  not  clear  which  index  i  gives  the  minimum,  except  that  the  likely 
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candidate  is  the  one  for  which  Y ^ ,  1  *  j  <  k#  is  a  maximum.  We  might  expect#  by 
analogy#  that  for  a  levels  of  practical  interest 


b  . ' ,  4  3  l  . 
g(x)  g(x) 


(4.8) 


5.  NUMERICAL  EXAMPLES 

For  k  =  2  and  selected  n,  y,  x,  a  -  .05  and  .10,  Table  I  gives  bg(J)<  ag(jj 

and  the  best  upper  bound,  u 

r  g(x) 


u  ~  =  min  u  (n.)  , 

9<X)  1<i<k  “  1 


(5.1) 


where 


“  “  Ju  (n  )(ni  ■  tt0i1'tt0il  +  ’> 
a  i 


(5.2) 


and  toi  are  defined  as  in  (4.6). 


TABLE  I. 

L-M  Approximations  and  ag(J) 


(n1 ,n2 ) 

(Y^Y.,) 

(x, ,x2 

(10,10) 

(1,2) 

(0,1) 

(10,10) 

(1,2) 

(0,1) 

( 10,10) 

(1,2) 

(1,1) 

(10,10) 

(1,2) 

(1,1) 

( 10,10) 

(1,2) 

(2,1) 

(10,10) 

(1,2) 

(2,1) 

(10,15) 

(2,3) 

(0,1) 

(10,15) 

(2,3) 

(0,1) 

(  10,15) 

(2,3) 

(1,1) 

(10,15) 

(2,3) 

(1,1) 

(10,15) 

(2,3) 

(2,1) 

( 10,15) 

(2,3) 

(2,0) 

a 

b 

g(x) 

*g(x> 

ug(x) 

.05 

.3670 

.3670 

.3670 

.10 

.4398 

.4398 

.4398 

.05 

.3045 

.3514 

.3670 

.10 

.3715 

.4227 

.4398 

.05 

.2484 

.3151 

.3670 

.10 

.3088 

.3825 

.4398 

.05 

.3695 

.3719 

.3742 

.10 

.4425 

.4446 

.4467 

.05 

.2554 

.3042 

.3670 

.10 

.3167 

.3705 

.4398 

.05 

.1712 

.1981 

.2431 

.10 

.2203 

.2513 

.3029 
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Mote  that  for  all  the  casea  In  T*ble  I,  b 

glx) 

computations  were  done  by  a  short  FORTRAN  program,  a 
the  author. 


is  a  lower  bound  for  agjjj.  The 
listing  of  which  can  be  obtained  from 


In  this 
components. 


6 .  CONCLUDING  REMARKS 

paper  we  have  extended  the  I/-M  method  to  series  systems  with  repeated 
More  work  is  needed  to  ascertain  the  region  of  validity  of  (4.6). 


i 

I 
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